Transport and Quantum Walk of Nonclassical Light in Coupled Waveguides 
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We study the transport and quantum walk of nonclassical light in an array of coupled waveguides 
which have novel properties like very low decoherence and thus making them ideal for storage of 
quantum information. We show how squeezing gets turned over from one waveguide to another. We 
further show how input nonclassical light can generate entanglement among different waveguides. 
Our results involve both first quantization due to array structure and second quantization due to the 
quantum nature of fields and can also be used to discuss the Talbot effect in the quantum regime. 
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Optical elements like beam splitters are known to be- 
have quite differently when it comes to single photons. A 
single photon according to Dirac can either be transmit- 
ted or reflected [lj . It can not be found simultaneously in 
both transmission or reflection. The behavior of a photon 
pair on a beam splitter is even more remarkable as shown 
by Hong, Ou and Mandel For a 50-50 beam splitter 
both photons are found either in transmission or in reflec- 
tion. This leads to entanglement of photons at the two 
output ports. The question of entanglement at the out- 
put of a beam splitter was investigated in very general 
terms [3, H|. Even more remarkable behavior has been 
shown experimentally d, 0] and theoretically Q . For 
other optical elements like phase shifters a quantum field 
with precise number of photons n undergoes a phase shift 
which is ncj) whereas a classical beam undergoes a phase 
shift <j). It is thus important to understand the behavior 
of single photons and more generally nonclassical light 
at other optical elements. Very recently coupled wave- 
guide systems, which are relatively easy to fabricate [9| 
and which are also relatively decoherence free, have been 
shown to be good candidates for continuous time random 
walks • The paper by Perets et al. [n| deals with clas- 
sical beams of light. In the light of what we have learnt 
with beam splitters and phase shifters and the fact that 
Feynman [111 ] used the term quantum walk to describe 
the behavior of quantum particles, it is worth examining 
how single photons and more generally nonclassical light 
would behave in coupled waveguide systems. This way 
we would be able to understand quantum walk by single 
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photons in coupled waveguides [12j, |13|, |14| . 

In this letter we consider the system of coupled wave- 
guides and report the propagation of single photons and 
nonclassical light. We discuss how nonclassical light from 
one waveguide gets transported to other waveguides. 
This is especially important in applications to quantum 
information science where one is very often interested in 
the storage and retrieval of a quantum state H| ■ We 
further report how nonclassical light at one input port 
can entangle different waveguides. We present analytical 
results for Heisenberg operators and wave functions for 
fields in different waveguides. For coupled waveguides, 



we show an analog of the well-known Hong-Ou-Mandel 
two-photon interference. We also report the amount of 
squeezing that can be produced in different waveguides. 

We consider an array of N single-mode waveguides, 
coupled through nearest-neighbor interaction. We will 
number the waveguides from 1 to N. The mode for the 
field in the j th waveguide is described by the annihilation 
(creation) operator a. (at). The operators and at for 
the coupled waveguide system obey the boson commuta- 
tion relation [a^aj] = 1. 

The Hamiltonian for the system can be written as, 

N N-l 

H = frg a j a j + hJ Yl («] a j+i + a j+i a j ■ ) • (!) 

3 = 1 3=1 

In the above equation, TV is the number of waveguides 
and the coupling parameter J represents the rate at 
which the photons are transferred to the neighboring 
waveguides. The Hamiltonian (1) can be diagonalized 
by using the normal co-ordinates given by 



N 

p=l 

N 
3=1 

where the function S(j,p) is defined as 
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S(j,P) 
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• sin 



N + l V7V + 1, 
This function satisfies the orthonormality relations 



(4) 



N 



^2S{n,p)S(m,p) = S nm , 



P =i 



N-l 



^2 (S(n,p)S(m,p + l) + S(n,p + l)S(m,pj) 



/ nix \ 



2 



These two relations lead to the diagonalization of (1) 



N 

H = h^ig + (3 P )blb p , 
p=l 



(G) 



Using (2)- (6) we obtain the Heisenberg operators for the 
field in each waveguide 



For input single photons in say waveguides i and j the 
probability of finding one photon in the waveguide k and 
one in the I th waveguide is \Ai^Aj t i+Ai^Aj^\ 2 ■ The two 
quantum amplitudes can obviously interfere. In particu- 
lar let us consider if the coupled waveguides can exhibit 
an analog of Hong-Ou-Mandcl interference. Consider two 
waveguides with one photon put in each, with a delay of 
T seconds. The wave function at time t can be related 
to the initial wave function using the evolution operator 
U(t): 



N 



Oi(t) = $>'(°) A . 



1 = 1 
N 



=^2exp[-i(g + P p )t]S(l,p) S{j,p) . (7) 
P =i 

The coefficients Aj i determine all the quantum prop- 
erties of light in the coupled waveguide system. For ex- 
ample if we start with light in the I th waveguide then 
\Aj i | 2 gives the propagation of light from the I th to the 
j th waveguide. The mean number of photons in the j 
waveguide would be 



(8) 



\m) = 



U(t-T)alu(T)a\\0,0) 
((l,0|L/t(T)a 2 4c/(T)|l,0))V2 



(10) 



In Eq. (fT0|) a} 2 corresponds to the addition of a pho- 
ton in the second waveguide at time T. Further the 
denominator in (10) arises as we have to insure the 
normalization of the wave function ci2U(T)a\\0, 0) at 
time T. Using U(t - T) = U(t)W(T) and the defi- 
nition of the Heisenberg operators a(t) = W{t)aU{t), 
the numerator in (10) simplifies to U{t)a\(T)a\\0,0) = 
4(T - t)a\(-t)U(t)\0, 0) = 4(T - t)a\(-t)\0, 0). 

Using the solution of Heisenberg equations in this nu- 
merator and using (10) we find that the probability of 
finding simultaneously one photon at each output at time 
tis(6 = Jt, 6 = JT) 
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FIG. 1: Behavior of the normalized intensity as a function of 
r (r = Jt/n). The input is injected into the first waveguide 
and the number of waveguides in the system is N = 6 

This result is similar to the result for classical propaga- 
tion. Clearly (8) is independent of g. In Fig. 1 we show 
the normalized intensity when the input is a single pho- 
ton state. Using the analytical solution (7) we can also 
build up the wave function of the field at time t. The 
form of the wave function is rather involved. For input 
Fock states the state at time t can be obtained in terms 
of the Aj/B given by Eq. : 

|m, n 2 , . . . , n N ) -> (E A "*L^ —\0, 0, . . . , 0) . (9) 



p(t,T) = |(1,1|^)| 2 , 

= cos 2 (20-0 o )/(l + sin 2 (fl o )) , (11) 

which shows the two photon interference dip at 2t — T = 
ir/2J depending on the length (proportional to t) of the 
waveguides and the delay time. For a given structure 
such a dip can be scanned by varying the delay time. 

For initial excitation in a single waveguide the number 
of photons in each waveguide does not depend on the 
quantum characteristics of the input field. We therefore 
examine next the squeezing and entanglement aspects of 
the radiation in different waveguides. We investigate the 
propagation of nonclassical light across the coupled wave- 
guides. We assume that squeezed light is coupled into the 
first waveguide. The input at the first waveguide is 

|C) = cxp(-^(e l Vi) 2 +e- 4 W))|0), 

(12) 

where r is the magnitude of squeezing and <fi is related to 
the orientation of the squeezing ellipse. In what follows 
we set g = 0. Its effect can always be incorporated by 
carrying out a simple rotation. 
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FIG. 2: Time evolution of the squeezing factor Sj(q) as a func- 
tion of t (r = Jt/iv). The input is in the first waveguide and 
the number of waveguides in the system is N = 5. The mag- 
nitude and the phase of the squeezing parameter are chosen 
as r — 0.7 and <f) — respectively. 



In order to study the transport of inter-waveguide 
squeezing, we introduce the quadrature operators for 
the j th waveguide given by qj = (a- + a])/V2 and 

Pj = (a - —a*-)/\/2i. We also define the squeezing factors 
Sj(q) = {Aqj) 2 ~ 1/2 and s (p) = {A Pj ) 2 ~ 1/2 . Thus 
squeezing occurs when one of these expressions becomes 
less than zero. Using Eq. (0 we get 




t 




0.5 



0.1 - 



-0.1 



-0.3 



I 

- 



1.4 



1.1 



li.S 



I 



IJ.5 



0.2 - 



-0.1 
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(13) 



FIG. 3: The top (bottom) part shows the variation of Sj(q) 
(sj(p)) as a function of r = Jt/n for j — 1, . . . , 5, smeared 
out in the j direction. The magnitude and the phase of the 
squeezing parameter are chosen as r = 0.7 and <j> = respec- 
tively. The number of waveguides in the system is N = 5. 



where — (+) sign is to be used for the quadrature q(p). 
In particular for a system of two waveguides, we have 



cos' 



'(Jt) sinh(r)( sinh(r) — cos(^) cosh(r)) , 



s 2 (p) cot 2 (Jt) 



(14) 



Clearly q-quadrature is initially squeezed if tanh(r) < 
cos(<p). Note that for two coupled waveguides we obtain 
complete transfer of squeezin g al beit from g-quadrature 
to p-quadrature for Jt — it/2 [191 ]. 

For a system of three waveguides, we get the following 
results 

si{q) = / cos 4 (-^=) , 
S3(q) = / sin 4 ( — \ 



(V2Jt) . 

f = sinh(r) ( sinh(r) — cos(^) cosh(r)) 



s 2 {p)= f -sin 2 {V2Jt) , 



(15) 



The inter-waveguide transfer of squeezing is gov- 
erned by the factors cos 4 (Jt/- v / 2), sin 4 ( Jt/y/2), and 
2sin 2 (Ji/\/2) cos 2 (Jt/v / 2)- Again, the g-quadrature is 
initially squeezed if tanh(r) < co$((f>). Also, for the case 
of three coupled waveguides we obtain complete trans- 
fer of squeezing from the first waveguide to the third 
waveguide for Jt = ir/\/2. In Fig. 2, we display the 



time evolution of the quadrature squeezing for the case 
of five waveguides. The negative values of s q (l), s q (3), 
and Sq(5) clearly demonstrate the squeezing in the q- 
quadrature. Fig. 3 shows the quadrature squeezing when 
the input is given to the middle waveguide. 

We next show that input of nonclassical light to one 
of the waveguides can produce pairwise entanglement be- 
tween the inter-waveguide modes. We use the well-known 
criterion for entanglement between two continuous vari- 



able systems [20|, |21|, |22(. As a measure of entangle- 
ment we examine the correlation between two waveguide 
modes. 



M(i, k) = 



\ a j a k 



<<M> 



(16) 



The negativity of M is a sufficient condition for entan- 
glement. For Gaussian states this is both necessary and 
sufficient. A calculation shows that the joint state of the 
coupled waveguides is Gaussian. We calculate M using 
the solution (8) and the state in (12). Before we show 
numerical results, we discuss analytical results for two 
and three waveguides. In particular, for a system of two 
waveguides, 

M(l,2) = -(sinh(2r)(tanhr-sin(2Jt) sin <f>)) (17) 

and thus entanglement occurs for sin(2 Jt) sin(</>) > 
tanh(r). 
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FIG. 4: The correlation function M(j, k) as a function of r 
(r = Jt/iv) for the case of six waveguides. The magnitude 
and the phase of the squeezing parameter are chosen as (a) 
r = 0.7 and (j> = 3n/2 ; (b) r = 0.6 and <j> — n respectively. 



For a system of three waveguides, we find the results 



M(l,2) = 



M(l,3) 



- cos I —7= 11(3 — cos 



((3-cos(V2Ji)) sinh 2 (r) 



2™ 

-\/2sm(0)smh(2r)sin(\/2Ji)) , (18) 
1 



(19) 



cos(^) sinh(2r) sin 2 (v / 2Jt) 
(3 + cos(2V2 Jt)) sinh 2 ( 



Clearly for <f> — tt, the first and third could be entangled. 
In Fig. 4, we show the time evolution of M(j, k) for the 
case of six waveguides. The negative values of M(l,j) 
(j = 1 to 6) in Figs. 4(a) and 4(b) clearly demonstrate 
the entanglement between the inter-waveguide modes. 

In conclusion, we discussed the continuous time quan- 
tum walk of quantum particles (photons) in a physical 
system consisting of N coupled waveguides. We showed 
that the quantum walk with a light source having strong 
quantum character can produce entanglement between 
different waveguides in the array. We also studied the 
transport of nonclassical light across the coupled wave- 
guides. We could investigate several other interesting 



nonclassical situations, for example effect of launching a 
distribution of entangled photon pairs into the array and 
we would then have an analogy of the Talbot effect in 
second quantized set-up. 
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